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Abstract 

We calculate the interaction-induced deformation of 
the Fermi surface in the two-dimensional Hubbard 
model within second order perturbation theory. Close 
to half-filling, interactions enhance anisotropics of the 
Fermi surface, but they never modify the topology of the 
Fermi surface in the weak coupling regime. 
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1 Introduction 

Since the discovery of high-temperature superconductiv- 
ity, the structure of low- lying single-particle excitations in 
two-dimensional interacting Fermi systems has attracted 
much interest ||l|, |2| . A key role is thereby played by the 
shape of the Fermi surface which determines the phase 
space for residual scattering processes and thus decay 
rates and Fermi liquid instabilities . In a recent Monte 
Carlo simulation of the two-dimensional Hubbard model 
Bulut, Scalapino and White Q have found that strong in- 
teractions may not only deform the Fermi surface of the 
non-interacting reference system, but may even lead to a 
different topology. While the Fermi surface of the non- 
interacting Hubbard model (with nearest-neighbor hop- 
ping) is always closed around the origin in k-space for 
densities n < 1 , the strongly interacting system exhibited 
a Fermi surface closed around (tt, tt) for densities close 
to but below half-filling. This result raises the interest- 
ing question whether such a behavior occurs only in a 
strong coupling regime. Zlatic, Schotte and Schliecker 1^ 
have recently argued that the Fermi surface topology can 
be changed by arbitrarily weak interactions in the limit 
n — > 1, since in that limit arbitrarily small deformations 
could lead to a topologically different shape. 

In the following wc will analyze the interaction-induced 



deformation of the Fermi surface in the two-dimensional 
Hubbard model within second order perturbation theory. 
We show that interactions enhance the anisotropy of the 
Fermi surface for densities close to half-filling, but they 
do not change the topology by deforming a (0, 0)-centered 
surface into a (tt, 7r)-centered one in the weak coupling 
regime. 

The deformation of the Fermi surface in two- 
dimensional Hubbard models has been studied already 
earlier by Schonhammer and Gunnarsson ||^ and by Os- 
sadnik [Q. The aim of the former authors was to show 
that Kohn-Sham Fermi surfaces are in general not exact. 
The behavior of deformations in the limit n — » 1, however, 
has not been addressed in these earlier works. 

2 Fermi Surface and Perturbation Expansion 

The Hubbard model Hamiltonian is given by 

H=-tY.Y. S^'- + ^ E "jT n,i (1) 

d-.j') 3 

where t is the nearest-neighbor hopping amplitude and U 
the on-site Coulomb repulsion. The operators cj^. (cjo-) 
create (annihilate) fermions with spin projection a on lat- 
tice site j and nja = Cj^Cja- 

The non-interacting band-structure for nearest- 
neighbor hopping on a square lattice is 



= —2t (cos kx + cos ky 



(2) 



This leads to a strictly convex Fermi surface centered 
around k — (0,0) for densities n < 1, to a diamond 
shaped surface at half-filling (n = 1), and a (tt, 7r)-centered 
Fermi surface for densities n > 1 (see Fig. |l]). 

The Fermi surface of the interacting system can be ob- 
tained from the one-particle Green function 



G(k,uj) 



(3) 
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Figure 2: Feynman diagrams contributing to the second 
order self-energy. 



Figure 1: The Fermi surface of the non- interacting system 
{JJ = 0) for various densities n. 

where I](k, w) is the self-energy. The chemical potential 
fj, controls the average particle density n. The equation 

Re G-\k, c^) = a - (fk - Ai) - Re S(k, Ck) = (4) 

determines the energy ^k = Ck ~ A* of coherent single- 
particle excitations (quasi-particles) in the interacting 
system . The Fermi surface is the set of those points in 
k-space where the excitation energy vanishes, i.e. ^kp = 0, 
and thus determined by 



Ai + Rel](kj.,0) = 



(5) 



According to the Luttinger theorem the volume en- 
closed by the Fermi surface is related to the particle den- 
sity by 

(^e(,-ek) (6) 

We now calculate the interaction induced deformation 
of the Fermi surface to second order in the coupling 
strength U . The particle density is kept fixed by a suit- 
able choice of an interaction-dependent chemical potential 
jjL = iJ,{n, U) — ^0 + S^j,, where = Mo(^) is the chemical 
potential corresponding to density n at J7 = 0. We denote 
the Fermi wave vectors of the non-interacting system by 
hpo and the deformation vectors hy Shp — hp — lipQ. 

To first order in U, the self-energy is a real constant: 
5]i(k, w) = Un/2. To keep the density fixed, the chemical 
potential has to be shifted accordingly by — Un/2, 
which cancels Ei completely in G. 

To second order in [/, two Feynman diagrams con- 
tribute to the self-energy of the Hubbard model (see Fig. 
|). The second one is a real constant which is completely 
cancelled by a corresponding shift of the chemical poten- 



tial.!^ '^he first diagram, however, leads to a k-dependent 
contribution and generates a Fermi surface deformation. 
Expanding Eqs. (^ and (^ with k^? = hpo + Slip and 
fJ- = fJ'O + 6 fi in powers of U and comparing the second 
order terms, one obtains the relation 



S'kp2 = S^l2 - S2(kFO,0) 



(7) 



for the deformation of the Fermi surface to second order 
in U, where the chemical potential shift is given by the 
Fermi surface average 



SfJ.2 = 



Jd^kS{el~fi) S2(k,0) 



(8) 



Note that S2 is real at zero frequency. Of course there 
are many ways to define a map k^g 1-^ hp between non- 
interacting and interacting Fermi surfaces. A natural 
choice is (5kf^ = (Jfe^nk^o, where nk^o is a unit vector 
normal to the non-interacting Fermi surface in hpQ. Eq. 
(^ then determines the modulus of the shift as 

(5^2 - S2(kF0,0) 



where v? 

KFO 



Skp2 = Skp2{kpo) = 



(9) 



j^^^ I is the Fermi velocity of the non- 
interacting system. (See also the equivalent expression 
(4) in §) 

The second order self-energy has not yet been evalu- 
ated by purely analytical means. Using the spectral rep- 
resentation of the non-interacting propagator Gq, one can 
write the imaginary part of the self-energy contribution 
associated with the first Feynman diagram in Fig. |2| as 



ImE2a(k,a;) 



sgn(u;) 



dv 



J i^ny JO 
Imno(q, ly) ImGo(k-q, cj-z/) (10) 

^We note that this diagram is "anomalous" in the sense of Kohn 
and Luttinger [pxt, i.e. it vanishes if the thermodynamic limit is 
taken strictly at T = 0, while it yields a finite contribution if the 
zero temperature limit is taken after the thermodynamic limit. 
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where the imaginary part of the non-interacting polarisa- 
tion bubble is given by 



0.001 



Imno(q, J^) = 



sgn(j/) f cPp 



TT J (27r)2 
ImGo(p, ijj') ImGo(p + q, uj' 



and ImGo (k, w) — — 7rsgn(cj) (5[tj — (e^ — /io)] ■ Performing 
the integrals over w', p^; and in (^0|) and ( |ll| ) analyti- 
cally (eliminating thus three (5- functions), one obtains the 
expression 

ImE2a(k,t^) = -sgn(a;) j dqy j dpy j dv 

e(/io-eO) -e(/io 



E 



Y(l-/2)(4sin2(g,/2)-.g2) 

with the functions 

/ = {v - uj - fio)/2t - cos{ky - qy) 
g = v/2t + cos{py + qy) - cos{py) 



(12) 



(13) 



The summation variables q^ and p|J in (|T^) are the roots 
of the set of equations 



cos{K-qx) = / 
2sin(pj,-|-(72;/2) sin((7j;/2) = g 



(14) 



The remaining three-fold integral is easily computed nu- 
merically (e.g. via a standard Monte-Carlo routine on a 
work-station). Note that the representation in (|l2|) is par- 
ticularly suitable for a high resolution of the low-energy 
limit w ^ 0, since the integration region shrinks with lo. 
The full self-energy function can be reconstructed from 
its imaginary part lmT,2a by a simple Hilbert transform 



i;2a(k, w) 



dio' 



|ImS]2a(k,C^0| 

bj ~ bj' + iO+sgn(a;) 



(15) 



The constant contribution Y,2b from the second diagram 
in Fig. |2| need not be calculated since it is completely 
cancelled by a corresponding shift 5^2b of the chemi- 
cal potential. Our numerical results for the second or- 
der self-energy agree with those published recently by 
Zlatic et al. j5| , who computed via a sequence of fast 
Fourier transforms. They also agree with earlier results 
by Schweitzer and CzychoU ||ll[ and by Ossadnik Q. 

3 Results for the Fermi Surface Deformation 

To discuss explicit results for the Fermi surface deforma- 
tion we introduce polar coordinates in k-space. Points on 
the Fermi surface of the non-interacting system are thus 
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Figure 3: Shifts of the Fermi surface dkp2/U^ as a func- 
tion of the angle cf) for various densities n. 



specified by an angle (j). Due to the discrete symmetries 
of the square lattice it is sufficient to consider angles be- 
tween and 45 degrees. In Fig. || we show results for 
Skp2/U'^ as a function of for various densities n. Here 
and in the following we set the hopping amplitude t ^ I. 
At low densities n < 0.6 weak interactions tend to com- 
pensate anisotropics of the non-interacting Fermi surface 
while at densities 0.7 < n < 1 anisotropics are further 
enhanced by interactions, as observed already earlier in 
unpublished work by Ossadnik Q]. The diamond shaped 
Fermi surface at half-filling is of course not affected at 
all by interactions, as a consequence of the particle-hole 
symmetry of the Hubbard model (|l|). More generally, 
particle-hole symmetry maps the Fermi surface at density 
n onto the surface for density 2 — n by a simple (tt, 7r)-shift 
in k-space (the diamond at half-filling is thereby mapped 
onto itself). 

Let us now clarify whether weak interactions can mod- 
ify the Fermi surface topology at densities close to half- 
filling. The first Fermi point that may reach the Brillouin 
zone boundary (and thus introduce a different topology) 
upon increasing U is obviously the one at = 0, since 
it is closer to the zone boundary than any other Fermi 
point already in the non-interacting system and in addi- 
tion SkF{(f>) is maximal for = (if n > 0.7). Quantita- 
tive information on the Fermi surface deformation in the 
"critical" regime n ^ 1 and (/> ^ is provided in Fig. |4[ 
where we have plotted 5kp2/U^ as a function of density 
for various small angles </>. Within numerical accuracy, 
Skp2{4') tends to for n ^ 1 for all 0, as expected from 
particle-hole symmetry and continuity. Finally, in Fig. |^ 
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we show the critical couphng strength Uc(n) that is re- 
quired to make k^^^ = + S]iF2 reach the BriUouin zone 
boundary (at the point (tt, 0)). Close to half-filling, Uc{n) 
behaves linearly as a function of density and extrapolates 
to a rather big finite value in the limit n — > 1. Hence, 
interactions do not modify the Fermi surface topology of 
the two-dimensional Hubbard model within the pertur- 
batively controlled weak coupling regime at any density. 

Zlatic et al. have calculated the critical coupling 
strength from second order perturbation theory for n — 
0.97 (only), where they obtained a smaller Uc (by a fac- 
tor of about 2) than we did, although our results for the 
self-energy agree. The discrepancy arises because these 
authors have determined the new Fermi surface by di- 
rectly solving Eq. (|^) with the second order self-energy, 
while we have expanded Skp to second order in U, which 
is the order we really control. In the small U limit both 
procedures yield the same shift to order J7^, but quanti- 
tative differences arise for finite U. The qualitative result 
that Uciji) remains finite in the limit n — )■ 1 is thereby 
not affected. 

In summary, we have calculated the interaction-in- 
duced deformation of the Fermi surface in the two- 
dimensional Hubbard model within second order pertur- 
bation theory. Close to half-filling, interactions enhance 
anisotropics of the Fermi surface, but they never modify 
the topology of the Fermi surface in the weak coupling 
regime. 
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